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Introduction 


Philip  Mandel  (1967)  has  defined  ship  maneuverinq  as 
"the  controlled  chanqe  or  retention  of  the  direction  of 
motion  of  a ship  and  its  speed  in  that  direction."  The 
study  of  ship  maneuvering  includes  the  problem  of  maintain- 
ing a fixed  heading  (course-keeping)  as  well  as  the  problem 
of  changing  the  ship's  heading  (steering). 

Traditionally,  these  problems  have  been  attacked  by 
assuming  that  all  of  the  hydrodynamic  forces  and  moments 
that  act  upon  the  hull  can  be  expressed  as  functions  of  the 
instantaneous  velocities  and  accelerations  of  surge,  heave, 
and  sway  and  the  instantaneous  annular  velocities  and  accel- 
erations of  roll,  pitch,  and  yaw.  These  assumed  hydrodynamic 
forces  and  moments  are  then  expanded  in  a Taylor  Series 
about  a uniform  forward  motion  and,  provided  that  the  devi- 
ations from  the  uniform  forward  motion  are  small,  only  the 
linear  terms  of  the  Taylor  expansions  are  retained.  This 
procedure  leads  to  a set  of  linearized  equations  of  motion 
that  provide  the  definitions  of  the  various  "stability 
derivatives"  as  well  as  the  basis  for  the  traditional  exper- 
imental techniques  used  in  their  evaluation.  Once  the 
stability  derivatives  for  a particular  ship  have  been  deter- 
mined, the  linearized  equations  of  motion  are  used  to  predict 
the  steerino  and  course-keeping  capabilities  of  the  ship. 

Since  this  traditional  approach  to  the  problem  assumes 
that  the  forces  and  moments  are  functions  only  of  the 
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instantaneous  values  of  velocities  and  accelerations,  any 
possibility  that  the  history  of  the  motion  might  affect  the 
present  situation  has  been  excluded.  A new  approach  to  the 
problem,  which  does  not  exclude  the  history  of  the  motion, 
was  introduced  by  Cummins  (1962).  Cummins'  approach  to  the 
problem,  which  was  improved  by  Ogilvie  (1964)  and  Lin  (1966), 
differs  from  the  traditional  approach  in  the  description  of 
the  hydrodynamic  forces  and  moments  acting  upon  the  ship. 

In  the  approach  of  Cummins,  Ogilvie,  and  Lin,  one 
expresses  the  hydrodynamic  forces  and  moments  as  pressure 
integrals  over  the  entire  wetted  surface.  Although  this 
approach  appears  to  be  much  sounder,  it  is  also  more  diffi- 
cult since  it  requires  a knowledge  of  the  pressure.  Even 
though  one  can  formulate  this  problem,  it  is  much  too  com- 
plex to  be  solved  at  present.  Fortunately,  the  authors 
mentioned  provide  us  with  a systematic  approximation  pro- 
cedure whereby  one  can  express  the  pressure  integral  in 
terms  of  the  same  velocities  and  accelerations  used  in  the 
traditional  method.  However,  the  form  of  this  expression 
is  not  the  same  as  that  given  by  the  traditional  method 
except  under  rather  special  circumstances.  The  primary 
difference  is  the  appearance  of  convolution  integrals  that 
allow  for  the  possibility  that  the  history  of  the  motion 
might  affect  the  present  situation,  a possibility  that 
cannot  be  handled  by  using  the  traditional  approach. 

There  is  reason  to  believe  that  this  memory  effect  is 
often  small  and  that  either  approach  will  yield  good  results 


in  the  prediction  of  many  standard  ship  maneuvers.  Recent 
papers  by  Fujino  and  Motora  (1975),  Nomoto  (1975),  and 
Fujino  (1975)  , support  this  opinion.  One  of  the  objectives 
of  this  project  is  to  determine  just  how  large  a role  is 
played  by  this  memory  effect. 

In  the  present  work,  after  a brief  description  of  the 
traditional  approach,  we  begin  with  the  linearized  expres- 
sions for  the  forces  and  moments  as  aiven  by  Lin.  From 
these  expressions,  we  develop  a method  of  experimentally 
determining  the  necessary  stability  coefficients.  The 
experimental  technique,  while  new  to  maneuvering  problems, 
is  well  known  in  other  areas  and  is  often  referred  to  as  an 
impulse- response  technique.  It  would  appear  that  this  new 
method  makes  a more  efficient  use  of  the  planar  motion 
mechanism. 

An  examination  of  the  special  case  of  regular  oscilla- 
tory motion  leads  to  relationships  between  the  stability 
derivatives  measured  in  traditional  experiments  and  the 
stability  coefficients  found  by  the  impulse-response  tech- 
nique. Both  types  of  experiments  have  been  performed  and 
the  results  are  presented  and  compared. 

Ship  maneuvers  that  correspond  to  a few  standard  rudder 
commands  have  been  calculated  by  usinc  the  traditional 
equations  of  motion  as  well  as  the  equations  recommended  by 
Lin.  A comparison  of  the  predictions  supports  the  opinion 
that  the  memory  effect  is  often  small. 


The  Traditional  Approach  to  the  Problem 


To  lay  the  groundwork  for  the  comparisons  which  follow, 
a short  description  of  the  traditional  approach  to  the 
determination  of  stability  derivatives  is  necessary.  A 
more  complete  description  can  be  found  in  Mandel  (1967). 


Equations  of  Motion 

The  coordinate  svstem  0 x v z (Figure  1)  is  fixed  in 

' O O'  o o 3 

space  with  z taken  vertically  downward,  the  (x  , y ) plane 
o ' o ' o 

coinciding  with  the  undisturbed  water  surface,  and  x taken 

o 

in  the  general  direction  of  the  motion  of  the  ship.  Then 
the  motion  of  the  ship  is  completely  described  by  the  posi- 
tion of  the  center  of  gravity  (x  (t) , v (t),  z (t))  and 

og  og  og 

the  heading  angle  ty(t).  In  such  a system,  we  can  write 
the  equations  of  motion  from  Newton's  Laws: 


mx  = X _ 
og  oT 


my  = Y m 
i og  oT 


I $ = N 
zg  og 


where  the  dots  indicate  derivatives  with  respect  to  time 


X _,  Y _ = total  force  in  x,y  direction 

oT  oT  1 

N = total  moment  about  vertical  axis  throuqh 

og 

the  center  of  gravity 

I = moment  of  inertia  about  vertical  axis 

zg 

through  the  center  of  gravity 


fa 


1 


Equations  (4)  are  nothing  more  than  the  Newtonian 
equations  of  motion  written  in  the  ship  coordinate  system. 
The  difficulties  with  which  one  is  faced  result  from  the 
inability  to  specify  X^Y™,  and  Nm. 

11  X 


Linearization  of  the  Equations  of  Motion 

The  forces  and  moments  of  the  right-hand  side  of  eqs. 
(4)  are  composed  of  several  terms  and  it  is  assumed  that 
we  can  separate  these  into  two  parts.  The  first  part  con- 
sists of  all  those  forces  and  moments  which  result  from 
perturbations  of  the  ship's  motion  about  its  mean.  The 
second  part  consists  of  all  other  external  forces  exerted 
upon  the  hull:  wind,  waves,  propeller,  rudder,  etc. 

We  can  write 


X 


= X + X, 


where 


T E 

X = forces  created  bv  small  motions 


X£  = external  forces  on  ship. 


It  is  now  assumed  that  X is  function  of  the  variables 
u,u  ,v,v,\p  $.  And  since  it  is  assumed  that  X results 
from  small  perturbations  about  the  mean  motion,  we  shall 
expand  X in  a Taylor  series  about  u = u^: 

X = X ( u , u , v , v , 'l> , l ) 

= X(uQ,  o,o,  o,o,o)  +Au|^+u|^+v-^|  + 


• V • -i  Y , Y 

+ v -r—  + i>  + — + (higher  order  terms) 


"v 


'<X 


where  Au  = u - u 
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Noting  that  X is  created  by  the  perturbations,  we  set 


X ( u ,0,0, 0,0,0)  = 0.  And  since  the  hull  is  symmetrical 
o 

about  the  (x,z)  plane,  X must  be  an  even  function  of 

• • ^x  9X  3X  3X  . . . . 

v,v,i ip,'<p.  Therefore  — i ~ ~ > -vr  must  all  equal  zero 


v V !> 


dX  • 

If  we  keep  onlv  first  order  terms,  then  X = iu  rr—  + u - — 

c i u * • 


where  all  partial  derivatives  are  to  be  evaluated  at 
( u = u , u = 0 , v = 0 , v = 0 , ’ii  = C , ; = 0 ). 


For  the  lateral  force  we  get 


„ 9Y  , • 9y 

Y = v •« — + v - — 

OV  d • 

V 


3Y  + 3Y 

) • 

V * 


where  Y(u  ,0,0, 0,0,0)  = 0 and  4^  , t-L  - 0 since  a change 

O ’U  - 

in  forward  motion  will  not  produce  a lateral  force  on  a 
svmmetrical  hull.  Similarly 


m - 3N  • dN  * 9N  , 3N 
N = v + v — + ^ : — - 

JV  j • <3  * 3 

V Ip  'll 


Introducing  the  yaw  rate  r = v , we  define  the 

stability  derivatives  Y = 4~  , N = 4^-  , Y = ~ , et; 
1 v 9v  v 9v  ' r 9r 


Since  the  derivatives  are  to  be  evaluated  at  (u  ,0,0,0, 0,0)  , 

o 

they  are  assumed  to  be  constants  which  may  depend  upon  the 
Froude  number  and  the  shape  of  the  hull  but  not  upon  the 
nature  of  the  motion,  so  long  as  the  motion  remains  small. 
Rewriting  eqs.  (4),  retaining  only  linear  terms,  we 


obtain 


(m  - X • ) u - X = X„ 
u u E 


r 


-8- 

(m-Y’)v-Y  v + (rax  - Y*)r  + (mu  - Y ) r = Y_ 
v v g r o r E 

(5) 

(mx  - N*)v  - N v + (I  - N*)r  + (mx  u - N )r  = N„ 
gv  v 2 r gorE 

Note  that  the  surge  equation  is  independent  of  sway  and  yaw 
rates,  i.e.  to  this  first-order  linear  approximation  the 
surge  equation  is  not  coupled  to  either  sway  or  yaw.  Since 
we  are  concerned  here  with  the  effects  of  sway  and  yaw,  we 
can  now  concentrate  on  the  last  two  equations  of  (5). 

Evaluation  of  the  Stability  Derivatives 

Although  it  is,  in  principle,  possible  to  evaluate  the 
stability  derivatives  from  theoretical  hydrodynamic  consid- 
erations, it  is  not  a simple  task  and  the  usual  procedure  is 
to  evaluate  them  experimentally,  The  technique  consists  of 
taking  a geometrically  similar  model  of  the  hull  and  forcing 
it  to  move  in  a known  trajectory  (so  that  v,v,r,r  are  all 
known)  and  measuring  the  externally  applied  forces  and 
moments.  Typically,  the  imposed  trajectory  is  a sinusoidal 
oscillation  about  the  mean  path,  first  oscillating  in  pure 
sway  and  then  in  pure  yaw.  We  refer  to  these  experiments 
as  regular-motion  tests. 

Case  A:  Pure  Sway 

Suppose  we  impose  the  motion 

v = v cos  ut,  r = 0,  u = u 
o o 

Then  the  measured  force  and  moment  will  also  be  sinusoidal 
functions  of  frequency  to  with  phase  angle  e 


o 


Y = Y cos  (wt  + £,)  = Y.^cos  tot  + Yout  sin  U)t 


in 


N = cos  (wt  + £.,)  = N-Lncos  wt  + N0ut  S"'‘n 


The  eauations  of  motion  for  sway  and  yaw  become 
(m  - Y • ) (-  v^wsin  wt'  - Yv^v0  cos  = Yin  COS  wt  + 


(6) 


Y sin  wt 
out 


(mx_  - N;.)  (-  vq  w sin  wt)  - Nv(vq  cos  wt)  = N.n  cos  wt  + 


g v 


N . sin  wt 
out 


Solving  for  the  stability  derivatives,  we  obtain 


Y . 

Y = - -HI 

v v 


out 

Y-  = — — + m 

v V U) 
o 

(7)  N. 

in 

N = " r; 

v 


N , 

N.  = + mx 

v v w g 


Case  B:  Pure  Yaw 

Let  r = r cos  wt,  v = o,  u 


= u , the  equations  of 


motion  become 


(mx  - Y*)  (-  r w sin  wt)  + (mu  - Y ) (rn  cos 
q r o o r u 


(8) 


Y cos  wt  + Y . sin  wt 
in  out 


(I  - N-  ) (~r0w  s:'-n  + u0  " cos  ut) 

N cos  wt  + N . sin  wt 
in  out 


10 


>< 


and 


(9) 


Y 

r 


Y* 

r 


N 


r 


N ' 

r 


Y 


m 


o 

Yout 

ro" 


mx 

1 


u 

o 


N . 
in 


N 


cut 

r 

o 


+ I 


mx 

g 


z 


u 

o 


Experimental  Technique 

The  planar  motion  mechanism  (PMM)  used  to  impose  the 
required  motions  of  pure  sway  and  pure  yaw  is  the  same  one 
used  by  Paulling  and  Wood  (1962)  and  a schematic  is  given 
in  Figure  2.  The  PMM  is  attached  to  the  towing  carriage 
and  two  rods  which  can  be  oscillated  independently  connect 
the  PMM  to  the  model.  Goodman  (1960)  provides  a more  de- 
tailed description  of  a PMM. 

The  model  is  attached  at  two  points,  forward  and  aft 
of  midships,  by  means  of  strain-gauge  dynamometers  used  to 
measure  the  lateral  forces.  A linear  potentiometer  connected 
to  the  forward  rod  measures  the  lateral  displacement. 

Since  the  vibration  of  the  carriage  produced  an  un- 
acceptable noise  level,  it  was  necessary  to  pass  all  three 
signals  through  matched  low-pass  filters.  Although  this 
resulted  in  a greatly  improved  signal- to-noise  ratio,  the 
problem  was  never  totally  eliminated.  All  three  signals 
were  simultaneously  recorded  on  a strip-chart  recorder. 
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To  impose  pure  sway  motion,  it  is  only  necessary  to 

set  the  PMM  so  that  both  rods  are  oscillating  in  phase  with 

each  other,  but  pure  yaw  is  not  so  simple.  To  produce  pure 

yaw,  it  is  necessary  that  the  forward  rod  lead  the  after  rod 

by  a phase  angle  a = 2 tan  ! — where  d is  one  half 

u 

o 

the  distance  between  the  rods.  Therefore,  any  change  of  for- 
ward speed  or  frequency  necessitates  a readjustment  of  the 
phase  angle. 

For  the  calculation  of  the  stability  derivatives, 
knowledge  of  the  amplitudes  of  the  forces  and  their  relative 
phase  angles  are  necessary.  These  quantities  can  be  read 
directly  from  the  strip-chart. 

Problems  with  the  Traditional  Method 

The  results  of  regular-motion  tests  are  dependent  upon 
the  frequency  of  the  oscillatory  motion,  i.e.  the  stability 
derivatives  are  not  constants,  as  was  presumed,  but  are 
functions  of  frequency.  This  frequency  dependence  has  been 
observed  by  Paulling  and  Wood  (1962) , van  Leeuven  (1964) , 
and  others  since  then.  For  the  calculation  of  ship 
maneuvers,  one  generally  uses  the  zero-frequency  value  of 
the  stability  derivatives,  but  for  the  study  of  motions  in 
rough  seas,  it  is  necessary  to  know  the  extent  to  which  the 
stability  derivatives  depend  upon  frequency.  But  in  any 
case,  it  becomes  necessary  to  perform  a large  number  of 
tests,  oscillating  the  model  at  different  frequencies. 
Paulling  and  Wood  performed  approximately  650  separate  test 
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runs  to  evaluate  the  stability  derivatives  of  a Mariner  class 
ship  at  four  Froude  numbers. 

As  one  attempts  to  evaluate  the  stability  derivatives 
at  lower  frequencies,  it  becomes  increasingly  difficult  to 
accurately  measure  the  forces  and  one  inevitably  reaches  a 
frequency  below  which  no  good  results  can  be  obtained.  Some 
researchers  have  been  further  limited  by  the  short  length 
of  their  towing  tank,  since  a lower  frequency  requires  a 
longer  test  section.  The  zero  value  of  the  stability  de- 
rivatives is  found  then  by  attempting  to  extrapolate  the 
values  found  at  hiaher  frequencies. 

Another  problem  which  limits  the  experimenter  is  the 
reflection  of  the  transverse  wave.  If  one  tests  at  too 
low  a forward  speed,  the  wave  created  by  the  oscillatory 
motion  can  reflect  off  the  sides  of  the  tank  and  interfere 
with  the  model.  It  has  been  our  experience  that  this  prob- 
lem makes  it  almost  impossible  to  get  good  results  for 
low-speed,  shallow-water  tests. 
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Transient-Motion  Approach  to  the  Problem 


A new  approach  to  the  problems  of  ship  maneuvering  is 
presented  and  this  leads  to  a new  method  of  determining 
the  stabilitv  derivatives. 


Equations  of  Motion 

We  shall  again  start  with  the  Newtonian  equations  of 
motion  written  in  the  ship  coordinate  system,  eqs.  (4). 
Then,  to  simplify  the  problem,  we  make  the  following 
assumptions:  (1)  wa  assume  that  the  ship  is  sailing  in 
smooth  water  where  the  only  disturbances  are  those  created 
by  the  ship;  (2)  we  assume  that  if  there  is  any  rolling, 
heaving,  or  pitching  of  the  ship,  the  interaction  with 
surge,  sway,  and  yaw  is  not  significant;  (3)  the  fluid 
is  assumed  to  be  inviscid  and  irrotational . As  a con- 
sequence of  these  assumptions,  the  water  can  act  upon  the 
ship  only  through  normal  pressure  and  we  can  write 

m ( u — V’ii  — x LI 2 * * S)  = / pn  dS  + X 

g J s x E s 

(10)  m(v  + u-g  + x ip)  = f pn_  dS  + Y 

o V I! 

I $ + mx  (v  + u']i)  = / p(xn  - yn  ) dS  + n„ 

nr  rr  * *■  t r ■*  v K 


p = pressure 


n , n = x,  y components  of  unit  normal  vector 
x y 


S = wetted  surface 


X„,  Y = external  forces  acting  on  the  hull 

ili 


L 
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The  main  difficulty  is  the  evaluation  of  the  pressure 
integrals  on  the  right-hand  side  of  eqs.  (10).  We  now 
require  some  systematic  approximation  procedure  . This 
problem  has  been  attacked  by  Cummins  (1962)  and  Ogilvie 
(1964),  and  later  Lin  (1966)  provided  us  with  a solid 
foundation  within  perturbation  theory.  Assuming  only  a 
continuous  velocity  field  and  small  deviations  from  a 
uniform  forward  motion  (and  small  disturbance  of  the  free 
surface  created  by  the  motions),  we  can  linearize  the 
equations  of  motion  as  follows: 


mu  = -U  u - 6 Au  - jAu(t-x)N  (x)dx  + X 


XX 


XX 


XX 


(11) 


m(v+u  i+  x ip)  = - u v-3  v - / v(t-x)  N (T)dx 

° g yy  yy  o yy 

oo 

- V p - f l (t-x)N  ( x ) dx+  Y 

y\p  y P n yP  E 


I if+mx  (v+u  ij)  = -u  ip-B  p - / (t-x)N  (x)  dt 
z g o pp  pip  i r pp 


f v ( t-x ) N (x ) dx  + N 

t py  e 


where,  consistent  with  the  small-motion  assumption,  thf 
non-linear  terms  on  the  left-hand  side  of  eqs.  (10)  have 
been  dropped. 

The  surge  equation  is  not  coupled  to  the  equations  of 

sway  and  yaw  and  will  not  be  discussed  any  further.  The 

external  force  Y_  and  moment  N„  contain  all  forces  and 
E c, 

moments  not  contained  in  the  pressure  integral.  These 
might  include  forces  caused  by  the  propeller,  rudder,  wind, 
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waves  and,  in  the  case  of  our  experiments,  the  planar  motion 
mechanism.  The  convolution  integrals  in  eqs.  (11)  represent 
the  effect  of  the  history  of  the  motion.  Defining  yaw  rate 
r = i p,  and  rearranging  terms,  the  sway  and  yaw  equations 
become 

OO 

(m+u  )v  +3  v +•  / v(t-x)N  (x)dx  + (mx  +u  )r 

yy  yy  n yy  g y'i> 


+ (mu  +6  ) r + r ( 1 

o y-i>  1 

o 


■)N  . (T)dT  = Y 

Y'l>  E 


(mx  +u  ) v +3  v + f v(t- 

g py  :p  y ^ 


N (_ ) dx  + ( I + u ) r 

vy  z yp 


+ (mx  u + 6 ) r + f r (t-x)  N ( x ) dx  = N 

g o Pp  ' 'pili  £. 


In  contrasting  eqs.  (12)  with  the  linearized  equations 
used  in  the  traditional  approach,  eqs.  (5),  the  major 
difference  appears  to  be  the  presence  of  the  convolution 
integrals.  This  means  simply  that  the  present  approach 
allows  for  the  possibility  that  the  history  of  the  motion 
affects  in  some  way  the  hydrodynamic  forces.  In  the 
traditional  approach,  the  forces  exerted  by  the  water  on 
the  hull  are  presumed  to  be  dependent  only  upon  the  instan- 
taneous values  of  the  motion  of  the  ship.  It  should  also 
be  pointed  out  that  we  arrived  at  eqs.  (12)  via  a system- 
atic approximation  scheme  with  its  foundation  in  perturbation 
theory  and  in  our  opinion  this  approach  is  sounder  than  that 
used  in  the  traditional  method. 

The  problem  now  consists  of  finding  a method  of 
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evaluating  the  constants  p . . , 6.  . and  functions  N . . (t ) . In 

il  il  il 

principle,  these  could  be  found  from  a theoretical  approach, 
but  we  shall  follow  the,  hopefully,  simpler  path  of  deter- 
mining them  experimentally.  For  this  we  shall  need  the 
Fourier  transform  of  the  equations  of  motion. 


Fourier  Transform  of  the  Equations  of  Motion 

The  Fourier  transform  of  a function  f(t)  can  be 
defined  as 

f(w)  = J"  f(t)  e"laJt  dt 

*'  — co 


where  a sufficient  condition  for  the  existence  of  the 
transform  is  that  f (t)  be  absolutely  integrable.  If  f(t)  =0 
for  t < o,  then  we  can  write 


f (<*>) 


/,0°  — ioit  r 00 

f(t)e  dt  = / f (t)  cos  ut  dt 

o o 

-i  /°°  f(t)  sin  cot  dt 
1 o 

f (oj)  - i f (co) 
c s 


The  Fourier  inversion  theorem  gives  us 

f(t)  = JitC*  f (,J))  eiu)t  dco 

^ co  r A a 

= — f ! f (co)  cos  cot  + f (co)  sin  cot 

o L c 5 


dco 


2 r°°  A 2 c 00  C 

— f (co)  cos  cut  dco  = — f (co)  sin  cot  do) 

J O C TT  1 0 S 


We  shall  also  use 


f (w)  = i co  f (co)  where  f (co)  = j _x f(t)  e dt 


-icot 


* 


17 


and  the  convolution  theorem: 

if  h(t)  = /_jo  g(t-T)  f ( t ) dx 

then  h(to)  = g(w)  • f ( uj ) 

If  we  take  the  Fourier  transform  of  the  linearized 
equations  of  motion,  eqs.  (12) , we  get  the  following  pair 
of  equations: 


i w (m+u  ) + £ 


M ( w ) I v ( w ) 
y v 


+ |icj(mx  + „ ) * mu  * * N (x)  r(w)  = Y (w) 

L g y,  o yv  yv  J e 


i u>  (nx  + „ ) 

g -y 


+ iw  ( I 


N U) I v(x) 


) + mx  u 
g o 


* N U)  : r (w)  = NJJ 

V V E 


We  now  define  the  following  "stability  coefficients' 
noting  that  they  are  all  functions  of  frequency  : 


= 3 + N (x) 


1 

yy 

yyc 

c 

( w ) = w ( m + 

H ) - 

N l 

V 

2 

yy 

yys 

c 

(uj)  -8  + 

N ( 

Lj) 

3 

Py 

'Pyc 

c 

(w)  = w (mx 

+ w , ) 

- N 

(w) 

4 

g 

<py 

Pys 

c 

(w)  = mu  + 

8 + 

N (w) 

5 

0 

yt 

yvc 

c 

(w)  = w (mx 

+ v . ) 

- N 

(w) 

s 

g 

y<P 

yvs 

T 
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c U) 

= 

mx  u + S 

+ 

N (to) 

7 

go  ipip 

[plpC 

C (U>) 

_ 

w (T  + u ) - 

N (<o) 

8 

z 

tj^S 

substituting  into 

eqs.  (11), 

we 

obtain 

A 

(C 

+ i 

C ) (v 

- i v ) + 

(C 

+ i C ) 

(r  - 

i r ) = 

Y - i Y 

i 

2 

c 

s 

5 

6 

c 

s 

c s 

(C 

+ i 

C ) (v 

- i v ) + 

(C 

+ i C ) 

(r  - 

i r ) = 

N - i N 

3 

4 

c 

s 

7 

8 

c 

s 

c s 

Separating 

real 

and  imaginary  parts , we 

obtain 

A 

A 

/A 

C 

(U3)  V 

(to)  + 

C 

(to)V  (to)  + 

C 

(u>)  r ( oj ) 

+ c 

(to)r  (to) 

= Y (to) 

l 

C 

2 

s 

5 

c 

6 

s 

c 

c 

A 

(to)  V 

(J  ~ 

c 

( GO  ) V ( 10 ) + 

C 

(to)r  (w) 

- c 

(to)r  (to) 

3 

< :* 

ii 

l 

s 

2 

c 

5 

s 

6 

c 

s 

(15) 

/A 

A 

C 

(w)  V 

( 10 ) + 

c 

(to)v  (to)  + 

C 

(<o)r  (to) 

+ c 

(co)r  (to) 

= N (to) 

3 

c 

4 

s 

c 

9 

s 

c 

C 

/A 

(u)  V 

(to)  - 

c 

(to)v  (to)  + 

c 

(to)r  (to) 

- C 

(to)r  (to) 

= N (to) 

3 

s 

4 

c 

7 

s 

8 

c 

s 

The  importance  of  eqs.  (15)  is  two-fold.  First,  as  we 
show  in  the  next  section,  these  equations  give  us  the  capa- 
bility of  evaluating  the  stability  coefficients  C ,C  , ,C  . 

: 2 s 

Secondly,  and  perhaps  more  importantly,  they  provide  a means 

of  evaluating  the  path  of  the  ship  given  the  external  forces 

and  moments.  Of  course,  this  can  also  be  accomplished  using 

the  equations  of  motion  (12)  if  the  constants  u. . 8.  . and 

ij  . i: 

functions  N. . ( t ) are  known.  But  the  evaluation  of  these 

constants  and  functions  requires  a knowledge  of  (or  an 
assumption  about)  the  behavior  of  the  stability  coefficients 
since,  for  example 


J 


as  u * 
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N (t)  = - / ( C (u>)  - i ) cos  u)t  du) 

yy  tt  ' o i yy 

by  the  inversion  theorem.  However,  for  ship  maneuvers, 

Y(w)  and  N(uj)  will  go  to  zero  for  w > some  to  ; then  if  we 

i 

know  C ,C  , ,C  for  o < to  < w , we  can  find  v(to)  and 

12  8 — 1 

r(co)  by  equations  (15)  and  ultimately,  v(t)  and  r(t)  by 
the  inversion  theorem. 

Evaluation  of  the  Stability  Coefficients 

The  experimental  evaluation  of  the  stability  coef- 
ficients is  accomplished  by  taking  a geometrically  similar 
model  of  the  hull  and  giving  it  an  impulsive  motion  such  that 
v(t)  and  r(t)  are  zero  before  t = 0 and  after  t = T.  Then 
the  infinite  Fourier  transform  can  be  replaced  by  the 
finite  Fourier  transform  for  0 £ t <_  T. 

Case  A:  Pure  Sway 

Suppose  r(t)  = o,  v(t)  = y(t) , where  y(t)  is  the 
lateral  displacement  of  the  model;  then  eqs.  (15)  become 

C (u))v  U)  + C (w)v  (w)  = Y (w) 

1 C 2 S C 

C (oj)v  (w)  - C (w)v  (w)  = Y (w) 

1 s 2 c s 

(16) 

C { w ) v ( u) ) +C  ( 0) ) v ( u> ) — N ( cu ) 

3 c w S C 

C (U))v  (oj)  - c (oj)  v (oi)  = N (w) 

3 S u C S 

If  y(t),  Y ( t ) , and  N(t)  are  measured  for  ° i.  t T , then 
their  Fourier  transforms  can  be  calculated.  Then  eqs.  (16) 
can  be  solved,  frequency  by  frequency,  as  four  simultaneous 
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linear  equations  in  the  four  unknowns  C ,C  ,C  ,C  . Note 

12  3*. 

that,  in  principle,  one  such  test  will  give  us  the  stability 
coefficients  for  all  frequencies. 

Case  B:  Combined  Sway  and  Yaw 

Once  the  coefficients  C ,C  ,C  ,C  are  known,  any  im- 

12  3 4 

pulsive  motion  which  combines  sway  and  yaw  will  enable  us 
to  find  the  remaining  coefficients  C ,C  ,C  ,C  . In  practice 

5 6 7 0 

we  set  the  two  supports  of  the  planar  motion  mechanism  to 

be  180°  out  of  phase.  Then  if  y (t)  is  the  position  of  the 

i 

forward  support  and  y (t)  the  position  of  the  after  support, 

we  obtain  the  following  results:  y (t)  = - y (t)  , 

2 1 

Uo  1 • 

and  v(t)  = - — y (t) , and  r(t)  = j y^ (t) . Measuring 
y (t)  , Y ( t ) , N { t ) we  can  calculate  r(to),  v(cj),  Y(to),  N(co). 

i 

Rewriting  (15)  with  all  known  coefficients  on  the  right- 
hand  side  , we  find 


C (to)r  (to)  + C (ui)r  (to) 
5 c 6 s 


Y (to) 
c 


C (to)  v (u)  - C (to)  v (to) 

1 C 2 s 


C (w)r  (to)  - C (<o)r  (to) 
5 s 6 c 


Y (to) 
s 


C (to)  v (to)  + C (co)v  (to) 
1 s 2 c 


C (<o)r  (to)  + C (to)r  (to)  = N (to) 
7 c 8 S C 


C (w)r  (to)  - C (co)r  (to)  = N (w) 
7 s a c s 


C ( to)v  (to)  - C (to)v  (to) 

3 c 4 S 


C ( to ) v (to)  + C (to)v  (ui) 

3 s 4 c 


These  equations  can  be  solved,  frequency  by  frequency, 
as  four  simultaneous  equations  in  the  four  unknowns  C ,C  , 

5 6 

C ,C  . Now  we  see  that,  in  principle,  we  need  only  one 

7 8 

sway  impulse  and  one  combined  sway  and  yaw  impulse  to  evalu- 
ate all  eight  coefficients  over  the  entire  frequency  range. 
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In  practice,  the  situation  is  not  that  simple  and  more  tests 
might  be  required  to  achieve  sufficient  accuracy.  This 
problem  will  be  discussed  in  a later  section. 

Experimental  Techniques 

The  model  was  again  attached  to  the  planar  motion 
mechanism,  the  only  alteration  being  the  disconnection  of 
the  electric  motor  so  that  manual  power  could  be  used.  The 
output  signals  from  the  strain-gauge  dynamometers  and  the 
linear  potentiometer  were  filtered  and  recorded  on  a 4 
channel  FM  tape  recorder  and  later  digitized  at  250  samples/ 
second . 

Operating  the  planar  motion  mechanism  manually,  the 
experimenter  provided  the  impetus  to  initiate  the  sway  or 
yaw  motion  and  the  mechanism  was  allowed  to  coast  to  a 
smooth  stop.  The  resulting  impulses  varied  considerably 
between  experimental  runs,  but  typically  had  a duration  of 
about  1 second  and  maximum  energy  at  1.25  Hz  and  a maximum 
lateral  displacement  of  1 inch.  A second  series  of  ex- 
periments was  run  with  a slower  pulse  of  4 seconds  duration 
and  maximum  energy  at  0.25  Hz.  This  was  the  slowest  pulse 
that  would  yield  forces  large  enough  to  be  accurately 
measured  with  our  equipment. 

Another  series  of  experiments  was  run,  during  which  we 
attempted  to  produce  pulses  which  approximate  a step 
function.  As  will  be  explained  later,  such  pulses  yield 
the  best  results  for  very  low  frequencies. 


r 
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Computer  programs  were  written  to  calculate  the  Fourier 
transforms  of  the  digitized  data  and  to  solve  eqs.  (16)  and 
^ (17)  for  the  stability  coef f icients . All  data  processing 

was  performed  on  the  University  of  California's  CDC  6400 
computer . 

The  Existence  of  the  Fourier  Transform 

A sufficient  condition  for  the  existence  of  the  Fourier 
transform  of  f(t)  is  that  f(t)  be  absolutely  integrable: 

/Zoo  ; f (t)  dt  < - • 

In  practice,  we  require  that  f(t)  be  zero  for  all  t < o 
and  return  to  zero  after  some  time  T > o.  For  case  A,  pure 
sway,  there  is  no  problem  since  any  pulse  of  finite  duration 
will  give  us  v(t),  Y(t),  and  N(t)  equal  to  zero  for  t < o 
and  t > T.  For  case  B,  however,  the  only  way  to  achieve 
this  is  to  have  the  centerline  of  the  model  coincident  with 
u^  both  before  and  after  the  pulse.  In  practice,  this  is 
difficult  to  achieve. 

If  we  allow  the  model  to  come  to  rest  with  some  non- 

A A /V 

zero  drift  angle,  then  v(co),  Y (ui)  / and  N ( tu ) will  all  be 
non-existent  since  v(t),  Y(t),  and  N(t)  will  reach  some 
non-zero  constant  value  for  all  t > T.  But  note  that 
v(t),  Y ( t) , and  N(t)  will  all  go  to  zero  for  t > T,  and 

^ A A 

therefore  v(co),  Y(w),  and  N(w)  all  exist.  If  we  take  the 
derivative  with  respect  to  time  of  the  equations  of  motion 
(12)  and  then  take  the  Fourier  transform  of  these  new 


T 
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equations,  we  find  that  we  can  still  use  eqs.  (17)  provided 

~ ~ ~ ~ i i 

that  we  replace  v(to),  r(ui),  Y ( co ) , N(aj)  with  v(uj),  r(co), 

A A 

Y(w)  , tl  (w)  . 

It  would  then  appear  that  we  must  differentiate  the 

recorded  data  before  taking  the  Fourier  transform.  To  see 

that  this  is  not  the  case,  consider  f(t)  = f for  all 

o 

t <_  o and  f(t)  = fT  for  t >_  T . Then  we  find 

i,  ^ r00  -iwt  c T i...  -iwt 

f(u>)  = I f(t)e  dt  = J f(t)e  dt 

J — CO  ‘ Q 


upon  integrating  by  parts  we  find 


f ( w)  = 


■iu)t 


f (t)  e j + iui  f f (t)  e 

° 

— - 


-lut 


dt 


rT  c/..-icot  , -icoT  , 

= iu)  J f(t)e  dt  + frr,  e - f 

O 1 o 


Since  we  can  then  calculate  f (to)  without  differentiating 
f(t)  and  since  we  can  rewrite  eqs.  (17)  in  terms  of  v(w), 
etc. , we  can  use  an  impulse  which  has  a non-existent 
Fourier  transform  without  additional  complexity  or  loss  of 
accuracy  provided  only  that  the  Fourier  transform  of  the 
derivative  exists.  Alternatively,  we  could  extend  the 
definition  of  the  Fourier  transform  to  include  such  a 
pulse  by  defining 


ei  ^ rT  c i x.  \ — iwt  , i if  — iwT  f > 

f u>  = j f(t  e dt  - - (f„e  - f ) 

J 0 !D  T o 

and  we  see  that 

f (w)  = iui  f (w)  . 


L 
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It  should  be  pointed  out  that  if  one  uses  this  extended 
definition  of  the  Fourier  transform,  it  will  be  necessary 
to  extend  the  Fourier  inversion  theorem  also.  We  have 


f(t)  - f (o ) = ft  f (x  ) dx 


dr  — [°°  £ ( cu ) cos  (jot  da) 
J 0 TT  1 0 C 


,t  2 r°°  • 

= / dx  — / f (aj)  sin 

J r\  tt  J n c? 


T\  ' 0 S 


got  da) 


Noting  that  f (oj)  = cof  ( co ) and  f (oj)  = - wf  (w) 

C S sc 


2 |.00  . t A 

f(t)  - f(o)  = — J dw  u)  f (u)cos  tor  di 

TT  1 0 1 n 13 


2_  r00 

TT  ■*  0 


(•CO  ^ 

J do ) / - ojf  ( a) ) sin  cot  dx 

' n ' O r* 


After  integrating  we  have  the  expended  Fourier  inversion 
theorem  : 

f(t)  - f (o)  = — f (oo)  sin  cot  dw 


- — / f (at)  (cos  cot  - 1)  do  . 
TT  0 c 


Effect  of  a Filter 

As  mentioned  earlier,  it  was  necessary  to  filter  the 
signals  in  order  to  improve  the  signal- to-noise  ratio.  If 
f(t)  represents  any  of  the  signals,  filtering  it  with  a 
linear  filter  is  equivalent  to  replacing  f(t)  with 


I(t)  = /°°  f(t-x)W(i)  dx 

o 


r 


1 


where  W(x)  depends  upon  the  characteristics  of  the  filter. 

> Then 

f(w)  = f (a) ) • W (uj) 

i 

An  examination  of  eqs.  (15)  shows  us  that  if  all  the  signals 
are  passed  through  identical  linear  filters,  the  W(w)  will 
cancel  out  and  no  accuracy  is  lost  due  to  the  filtering. 

A Difficulty  in  Transient  Experiments 

As  mentioned  earlier,  since  a finite  pulse  has  com- 
ponents at  all  frequencies,  it  is  theoretically  possible 
to  run  one  sway  test  and  one  combined  sway  and  yaw  test 
and,  from  this  data,  solve  for  the  stability  coefficients 
over  the  entire  range  of  frequencies  o <_  o < °°  . But, 
since  we  are  passing  the  signals  through  a low-pass  filter, 
we  cannot  reasonably  expect  to  obtain  accurate  results  for 
frequencies  above  the  cut-off  frequency  of  the  filter, 

5 Hz  in  our  case.  This  is  not  a severe  limitation  however, 
since  5 Hz  is  a considerably  higher  frequency  than  one  needs 
for  almost  any  application. 

Unfortunately  we  face  a more  serious  problem.  When  one 
solves  eqs.  (16)  and  (17)  for  the  stability  coefficients, 

one  finds  expressions  for  C ,C  , ,C  Which  alwavs  contain 

12  8 

A A A 

a terra  in  the  denominator  such  as  (v  2 + v 2 ) or  (r  2 + r 2 ) . 

c s c s 


For  example: 
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Suppose  we  approximate  the  pulse  (for  case  A:  pure 

1 2 TT 

s' way)  by  v(t)  = =-(l  - cos  t),  o < t < T = — . We 

2 O — GO 

O 

shall  refer  to  such  a pulse  as  a "full  pulse".  Then 


v(u>)  = y(u>)  = 

1 

0 

. 

2tt  w 

2 

2 

OJ  " “ GJ 

0 

v z-  S X Xl  (>  — — . 

0) 

0 

• i sin 

GJ 

0 

Notice  that  v (u) 
c 

has 

zeros  at 

Cl ) ■—  O / 2 OJ  / 
0 

3<j  , 

0 

and 

v (to)  has  zeros  at 
s 

0)  = 

o , iu  , 

2 o 

, 2w  , 

2 o o 

. Then 

the 

denominator  in  eq.  (18) , (v  2 + v 2 ) , has  double  zeros  at 

c s 

cu  = o,  2w  , 3w  , . There  will  exist  a singularity 

o o 

at  these  points  unless  the  numerator  has  matching  zeros 

to  cancel  it  out.  The  terms  in  v in  the  numerator  orovide 

s 

only  simple  zeros,  but  one  assumes  that  the  transforms  of 
the  measured  forces  will  supply  the  additional  matching 

A 

ones.  In  practice,  this  cannot  be  realized  since  Y(w) 
contains  the  transform  of  the  signal  plus  the  transform  of 
the  noise,  and  there  is  no  reason  to  suppose  that  the  trans- 
form of  the  noise  goes  to  zero  at  these  frequencies. 

However,  even  if  this  were  so,  one  is  still  in  the  position 
of  dividing  two  very  small  quantities  at  and  in  the  neighbor- 
hood of  the  zeros  and  consequently  one  is  very  vulnerable 
to  small  errors  in  measurement,  which  may  become  very  large 
relative  to  the  quantities  measured. 

With  such  a pulse,  one  cannot  avoid  the  problem  at 
ij>  = o,  but  one  can  choose  a large  enough  u>Q  so  that  the 
remaining  zeros  are  outside  the  range  of  interest.  For 


i 
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example,  if  the  duration  of  the  pulse  T = 1 second,  then 
uj  ^ = 2*t  and  we  can  expect  reasonable  results  for  o < w < 4 tt  , 

but  if  T = 4 seconds  we  can  expect  reasonable  results  only 

* * 

for  o < u)  < tt  . Figure  3 gives  the  value  of  (v  2 + v 2)  for 

c s 

these  two  pulses. 

One  way  to  avoid  the  problem  at  go  = o might  be  to 

select  a pulse  which  does  not  return  to  zero,  such  as  y(t)  = 

1 TT  TT 

= ^-(1  - cos  a)  t)  , o < t < — where  v(t)  = 1 for  t > — . 

2 0 — — CO  GO 

0 0 

We  shall  refer  to  such  a pulse  as  a "step  pulse".  Then 

A _ UJ 

/\  */\  1 . , TT  CJ  • . TT  UJ  . 

V (lo)  = y (uj ) = - — (1  -f-  cos  — - l sm  — ) 

111  UJ  d) 


Notice  that  v 

c 

(w) 

has  zeros  at 

CJ  = OJ 

0 

t 3u  f 3 uj  f ~~ 

0 0 

— and 

that  v (w)  has 
s 

zeros  at  j = o, 

2uj  , 3 

0 

uj  , .Then 

0 

the 

denominator  in 

eq. 

(13) , (vc2  + 

V}' 

has  double 

zeros 

at  u>  = 3w  , 5 go  , 7go  , . Therefore,  such  a "step  pulse" 

o o o 

should  provide  good  results  for  o _<  io  < 3oj  . The  denomin- 
ator for  a step  pulse  with  coQ  = 2^  is  shown  in  Figure  3. 

A comparison  of  the  full  pulse  and  the  step  pulse  as 
used  in  our  experiments  is  shown  in  Figure  3a. 

There  is  another  way  to  avoid  this  difficulty,  although 
it  has  the  disadvantage  of  requiring  more  experiments. 
Suppose  we  perform  the  same  maneuver  several  times  with 
slightly  different  values  of  w , sav  w < w < to  . Then 

1 2 3 

we  can  sum  the  results  of  the  individual  runs  to  form 


- 

a r> 

2 u 

~ 

/V 

A 

A 

V (uj) 

= V 

+ 

V 

+ 

V 

1 

2 

3 

A 

A 

A 

A 

Y(u>) 

= Y 

+ 

Y 

+ 

Y 

l 

2 

3 

A 

A 

-A 

N (w) 

= N 

+ 

N 

+ 

N 

1 2 3 

If  the  three  runs  have  only  slightly  different  values  of 
w0 , one  can  show  that  (vc‘  + vs2)  is  not  likely  to  have 
any  zeros  near  2u)q  . This  procedure  has  been  followed  and 
Figure  4 compares  results  obtained  from  individual  runs 
and  the  result  of  the  combined  runs. 
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Relationship  Between  Traditional  Method  and  Transient  Method 


Note  that  the  equations  of  motion  used  in  the  impulse 
test  assume  nothing  about  the  motion  other  than  the  require- 
ment that  the  motion  be  small  perturbations  about  a uniform 
motion.  It  is  of  interest  to  examine  the  case  used  in  the 
traditional  method,  i.e.  regular-oscillatory  motion  about  a 

uniform  forward  speed.  Let  v = v coswt,  u = u , r = 0 and 

o o 

substitute  into  eqs.  (12): 


(m+u  ) ( - v wsinwt)  + 8 (v  coswt) 
yy  o yy  o 


+ / v cosw(t-x)N  (x)dt  = Y 
n o yy  E 


(mx  +u  , ) (—v  wsiniot)  + B . (v  cosuit) 
g ipy  o {iy  o 

oo 

+ f v cosw(t-x)N  (x)dx  = N 
o o :l>  y E 


Examining  the  convolution  integral,  we  find 


I v cosuj(t-")N  (x)dx  = / v coswtcosuix  N (x)dx 


+ f v sinwt  sinuixN  (x)dx 

o ° yy 


= v cosut  N ( uj ) tv  sinwt  N (w) 
o yyc  o yys 


Using  this  relationship  and  separating  the  force  and  moment 
into  their  in  phase  and  out  of  phase  components,  we  obtain 
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[oj(m+y  )-N  (to)  ] (-v  sinojt)  + [B  +N  (u)]v  cosut 
yy  yys  o yy  yvc  o 


/in.  = Y.  coswt  + Y sirn^t 

(19)  in  out 


[o)(mx  +y  ) -N  (u)  ] (-v  sinjtJ  + 18,,  +N , (ui)]v  cosot 


g ipy  i(iys 


ipy  ipyc 


= N.  cosot  + N sincjt 
in  out 


Comparison  with  eqs.  (6)  gives  the  following  relationships : 


Y = -8  - N (o) 

v yy  vyc 


v _ , -iN  (o) 

Y*  - -y  +o)  vys 
v yy 


(20) 


..  _ a ~ N (<*>) 

N = ~p , vyc 

v iv  1 


N * = - U + ‘ N ( o ) 

v yy  -ys 


and  a similar  examination  of  oscillatory  yaw  motion  yields 
Y = 


■8  - N (w) 

yy  y-yc 


Y-  = 
r 


(21) 


N = 
r 


N * = 
r 


'"‘yll) 

+ 0)  N 

yips 

_6yy 

3 

U 

■3- 

-B- 

< 

i 

- u 

+ u)~  1 N , 

y^s 

Finally,  a comparison  with  eqs.  (14)  yields 


C (id)  = - Y 
l ” 


C ( u> ) = U)  (m-Y  • ) 

•>  V 


C (u)  = -N 
3 v 


J 
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c 

(w) 

= to  (mx  - N • 

4 

g v 

c 

( (Jj ) 

= mu  - Y 

5 

o r 

c 

(w) 

= io  (mx  - Y • 

6 

g r 

c 

(oo) 

= mx  u - N 

7 

go  r 

c 

(to) 

= to  (I  - N-  ) 

8 

z r 

If,  as  presumed  by  the  traditional  method,  the 

stability  derivatives  are  constants,  then  N (to)  must  be 

yyc 

A 

a constant  and  N (to)  must  be  zero  everywhere.  This  means 
yys 

that  N (x)  can  be  written  as  a delta  function 


N (t)  = N 6 ( i ) where  N = constant 

yy  yy  yy 


Then  the  convolution  integrals  appearing  in  the  equations 
of  motion  can  be  written  as 

CO 

/ v(t-x)N  5(x)dx  = N v ( t ) 

o yy  yy 

and  any  dependence  upon  the  history  of  the  motion  is  lost 
and  indeed  the  equations  of  motion  take  on  a form  which  is 
identical  to  that  used  in  the  traditional  approach.  There 
is,  then  an  equivalence  between  the  dependency  of  the 
stability  derivatives  upon  the  frequency  of  oscillation  and 
the  dependency  of  the  instantaneous  forces  and  moments  upon 
the  history  of  motion.  Therefore,  the  fact  that  previous 
studies  have  shown  that  the  stability  derivatives  are 
frequency  dependent,  forces  one  to  conclude  that  the 


traditional  equations  of  motion  are  not  adequate  to  describe 
all  situations  and  that  the  convolution  integrals  should  be 
included  in  the  equations. 

Equations  (22)  give  us  a means  of  comparing  the  results 
of  the  two  experimental  techniques,  i.e.  regular-motion 
tests  vs.  impulse  tests.  Therefore,  rather  than  present 
the  results  in  terms  of  the  stability  coefficients 


C , C , , C , we  chose  to  present  everything  in  terms  of 

1 2 8 


the  more  familiar  stability  derivatives  Y , Y*,  , N*  . 

v v r 


1 


Low- Frequency  Behavior 

Noting  that  the  Fourier  cosine  transform  is  always  an 
even  function  of  frequency  and  that  the  sine  transform  is 
odd,  an  inspection  of  eqs.  (20)  and  (21)  leads  us  to  con- 
clude that  all  the  stability  derivatives  must  be  even 
functions.  Therefore,  if  we  express  the  stability  deriva- 
tives as  a Taylor  expansion  about  cu  = o,  we  have 

82Y 

Y = Y (uj=o ) + i + 

v 2 , 2 

do) 

It  is  now  apparent  that,  when  one  attempts  to  extrapolate 

regular-motion  test  results  to  qj  = o,  one  may  assume  that 

the  stability  derivatives  approach  a constant  value  with 

zero  slope.  Furthermore,  there  must  exist  some  range  of 

frequencies  o w < e over  which  the  approximation 

Y = Y (w=o)  is  usable.  Obviouslv,  if  the  forces  and 
v v 

moments  applied  to  the  hull  have  frequency  components  which 
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are  primarily  within  this  range,  then  there  should  be  little 
error  in  using  the  traditional  aoproach. 

Almost  Steady  Motion 

Since  it  has  been  shown  that  the  traditional  approach 
to  maneuvering  problems  will  yield  reasonable  results  for 
many  standard  ship  maneuvers,  it  is  of  interest  to  examine 
the  conditions  under  which  the  traditional  equations  of 
motion  (5)  become  a good  approximation  to  the  preferred 
equations  of  motior  (12).  For  the  sake  of  simplicity,  con- 
sider the  sway  equation  for  the  case  where  r = o.  Then 
eq.  (5)  becomes 

(23)  (m-Y-)v-Y  v=Y 

v v E 

where  Y*  and  Y are  to  be  evaluated  at  to  = o.  Note  that 
v v 

there  is  nothing  in  the  derivation  of  eqs.  (5)  which  allows 
one  to  assume  that  the  zero-frequency  value  of  the  stability 
derivatives  should  be  used.  However,  when  examining  ship 
maneuvers,  one  is  dealing  with  very  slow  motions  which 
suggest  a similarity  to  very  low  frequency  regular-motion 
tests . 

Taking  the  zero-frequency  limit  of  eqs.  (20) , we  see 

that 


(24) 


Y = -0 

v yy 


r 00 


M 

yy 


(x ) dt 


Y.  = -y  + r TN  (T)dr 
v Myy  ]o  yy 


3 4 


Substituting  eqs.  (24)  into  the  traditional  sway  equation 
(23) , we  find 


(25) 


(m  + u - / xN  (x)dx)v  * (S  + f°°  N (x)dx)v 
yy  o yy  yy  ; 0 yy 


= Y. 


The  transient-motion  swav  eauation  is 


(26) 


(m  + u )v  +3  v + v(t-x)N  (r)dx  = Y 
YY  yy  o yy 


We  now  ask  ourselves,  "Under  what  conditions  will  the 
solution  of  eq.  (25)  be  a good  approximation  to  the  solution 
of  eq.  (26) ?" 

Let  v (t)  be  the  solution  of  eq.  (26)  and  v (t)  be  the 
1 2 

solution  of  eq.  (25),  given  the  same  initial  conditions  and 
forcing  functions  for  each.  Subtracting  eq . (25)  from 

eq.  (26)  and  rearranging  terms,  we  find 


(27) 


(m  + u - f IN  (x)d-r)  (v  - v ) 

yy  0 yy  i 2 


+ (S  + / N (x)dx) (v  - v ) = g ( t ) 

yy  0 yy  1 2 


where 

(29) 


g (t) 


f0  ir(v  (t) 

0 L:  1 


v ( t-x  ))  - v ( t)  1 tN  dx 

1 1 yy 


If  we  assume  that  v = v at  t = o and  if  we  define 

i 2 

the  error  E = v - v , then  the  solution  to  eq.  (27)  is  of 

1 2 


the  form 


Apparently,  the  smaller  one  chooses  e (smaller  acceptable 
error) , the  larger  the  value  of  t will  become.  Let  us 
examine  the  result  of  placing  the  following  restrictions 
upon  the  acceleration: 

(32)  | v(t)  j < j ^ for  all  t 

(33)  | v ( t ) - v(t-t)  | < for  T 1 fcc 

By  the  mean-value  theorem 

v (t)  - v (t  - t)  = t v (t  - at) 

i i i 

where  0 < a(t,r)  < 1,  so  that  the  definition  of  g(t), 
eq.  (28) , can  be  rewritten  as 

•—  *»  — 1 

(34)  g(t)  = f | v ( t-cx  t ) - v (t)  t N dT 

O 1_  1 1 _J  n 
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1 


L 


Making  use  of  eq. (33) , we  see  that 

l»,l  < »yyldT  < 7£  K Nyy  | dT 

and  using  eqs.  (31)  and  (32) , we  can  see  that 


i V r°°  i , , 1 V f c 

gJ  < T ft  lTNyy'dT  = I e T f 


T N ! dx 

yy 


or 


1 V z*00  I 1 

q < — £ — I I t N d i 
1 2 2 T'  o 1 yy 


Therefore 


or 


g(t) | < £ J /"|t  Myy!dx 


g ( t)  | < eV  ( 3 + jT  N df) 

yy  0 yy 


from  the  definition  of  T.  Substituting  the  maximum  value  of 
g(t) ( into  eq.  (20)  we  reach  the  following  result: 

‘ E | • £ V (l-e7t)  < £V  . 

Therefore,  if  the  maximum  acceptable  error  is  to  be  eV, 
and  if  the  acceleration  meets  the  requirements  of  eqs.  (32) 
and  (33),  then  one  may  use  the  traditional  equations  of 
motion  with  the  coefficients  evaluated  at  jj  = 0. 

A different  approach  to  this  problem  can  be  found  in 
Wehausen  et  al.  (1976). 
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A large  number  of  experiments  have  been  performed  at 
the  University  of  California  in  an  attempt  to  perfect  the 
impulse-response  technique.  This  section  outlines  the 
various  attempts  which  led  to  the  currently  favored  method. 

In  all  cases,  the  experiments  were  performed  at  the 
University's  Richmond  Field  Station.  The  towing-tank  is 
approximately  200  feet  in  length,  8 feet  wide,  and  6 feet 
deep  (the  water  level  was  maintained  at  the  maximum  depth 
throughout  the  experiments).  The  planar-motion  mechanism 
which  was  used  is  the  same  one  used  by  Pauli. r.  : ■. n : Woe.: 

(1962).  The  model  that  was  used  is  a ..  :ht-v.  . * 

model  of  a high-speed  ship  (DE  type  ar.  : . • * : e 

tests  were  performed  using  the  model  w : t : .t  •••>;  •• . r 

rudder.  It  should  be  pointed  out  t:.  it  t.-..  : : . . 

propeller  and  rudder  in  no  way  affect  t..»>  -ta. 

technique  and  in  fact  another  researcher  ner»  . : . as 

Loeser)  has  performed  impulse  tests  usm  ; a Mariner  model 
equipped  with  propeller  and  rudder  [see  Wehausen  et  al. 
(1976)].  The  dimensions  of  the  model  are  as  follows: 

L = 5.0  feet 

B = 0.585  feet 

T = 0 . 19  feet 

C = 0.492 

B 

M 
I 


2 


0.239  slugs 
0.468  slug  • ft: 
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As  mentioned  earlier,  eqs.  (22)  give  us  the  ability 
to  present  the  results  of  both  regular-motion  tests  and  im- 
pulse tests  in  terms  of  either  the  traditional  stability 

derivatives,  Y , Y • , N , or  the  stability  coefficients 

v v v 

C , C , C , . . . In  order  that  the  present  results  might  be 

1 2 3 

more  easily  compared  with  the  work  of  other  researchers,  we 
present  all  results  in  terms  of  the  traditional  stability 
derivatives.  The  stability  derivatives  are  made  dimension- 
less with  i-P  , L,  and  u following  the  "prime  system"  used 

2 o 

by  Mandel  (1967).  Two  dimensionless  forms  of  the  frequency 
are  used: 

do  U 

_ o 

g 

and  ' _ wL 

u 

o 

We  note  that  t = F 2 w ' . 

n 

Regular-Motion  Tests 

In  order  that  we  would  have  data  with  which  to  compare 
the  results  of  the  impulse- tests , it  was  necessary  to  per- 
form a number  of  regular-motion  tests.  These  experiments 
were  performed  by  Tomas  Frank  (1974)  using  traditional 
planar-motion  mechanism  techniques.  Since  each  regular-motion 
experiment  yields  the  value  of  the  added  mass  and  damping 
coefficients  at  one  particular  frequency,  the  results  of 
these  experiments  appear  as  individual  data  points  and  no 
attempt  at  curve- fitting  has  been  made. 

Note  that  there  is  a range  of  low  frequencies  (see 


1 
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Figures  5-12)  in  which  no  results  are  given.  This  is  an 
inherent  problem  of  regular-motion  testing.  As  mentioned 
earlier,  one  inevitably  reaches  some  frequency  below  which 
accurate  measurements  are  impossible. 

Due  to  limitations  on  the  accuracy  of  the  measurements, 
the  results  of  these  tests  are  likely  to  contain  errors  on 
the  order  of  10-15  per  cent  and  therefore,  in  the  compari- 
sons which  follow,  the  differences  between  the  results  of 
the  two  experimental  techniques  should  not  be  regarded  as 
a measure  of  the  inaccuracy  of  the  impulse-test  procedure. 


Full-Pulse  Impulse  Tests 

The  first  series  of  experiments  from  which  we  received 
reasonable  results  employed  a full-pulse,  as  described 
earlier  (Figure  3a) , with  a duration  of  approximately  one 
second  and  a lateral  displacement  of  one  inch.  Such  a pulse 
has  its  peak  energy  at  about  1.25  Hz  and  will  yield 
reasonable  results  for  some  range  of  frequencies  centered 
about  this  point.  In  a previous  section,  it  was  explained 
that  such  a pulse  will  lead  to  results  which  are  singular 
at  u = 0,  and  indeed  this  problem  was  encountered. 

In  an  attempt  to  obtain  better  results  at  lower  fre- 
quencies, a second  series  of  experiments  was  run.  This 
time,  we  used  the  longest-duration  pulse  for  which  we  could 


still  measure  the  forces  accurately  with  our  equipment. 
These  pulses  averaged  four  seconds  in  duration,  peak  energy 
at  0.25  Hz,  and  one  inch  lateral  displacement. 
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Figures  5 through  12  show  the  results  of  these  two 
series  of  experiments,  as  well  as  the  results  of  the  regular- 
motion  tests,  for  = 0.30.  The  graphs  show  that,  though 
the  longer-duration  pulses  did  yield  slightly  better  results 
for  low  frequencies,  the  improvement  was  limited  to  a dis- 
appointingly narrow  range  of  frequencies.  It  also  becomes 
evident  that,  if  one  desires  information  about  the  zero- 
frequency  limit,  a different  sort  of  pulse  is  required. 

The  Step-Pulse  Impulse  Tests 

As  was  pointed  out  in  an  earlier  section,  a pulse  which 

approximates  a step- function  does  not  have  the  problem  of 

singularities  at  o = 0.  Therefore,  another  series  of 

experiments  was  run  using  the  step-pulse  (see  Figure  3a) . 

The  results  of  these  experiments  (with  an  inch  displacement) 

are  presented  in  Figures  13  through  20  for  F^  = 0.30  and 

Figures  21  through  28  for  F = 0.20. 

n 

These  graphs  indicate  that  all  of  the  damping  coeffi- 
cients are  well  behaved  at  w = 0.  For  the  case  of  pure 
sway,  the  added  masses  are  similarly  well  behaved.  However, 

the  two  added-mass  terms  Y*  and  N*,  which  are  calculated 

r r 

from  the  case  of  combined  sway  and  yaw,  still  "blow-up"  for 
u = 0.  Therefore,  a further  examination  of  this  case  appears 
necessary. 

The  Zero-Frequency  Correction 

In  order  to  understand  the  behavior  of  Y*  and  N • at 

r r 


"T 


zero  frequency,  it  is  necessary  to  return  to  eqs.  (17)  from 

which  we  calculated  C and  C , the  corresponding  stability 

6 8 

coefficients.  For  the  sake  of  simplicity,  only  C will  be 

6 

examined  here,  since  the  examination  of  C follows  a similar 

8 

path. 


Since  the  step-pulse  leads  to  a non-zero  force  before 
and  after  the  pulse,  we  shall  deal  with  the  Fourier  trans- 
form of  the  derivative  of  the  force,  which  exists  in  the 
conventional  sense.  Rewriting  the  first  pair  of  eqs.  (17), 
we  obtain 

C (w)r  (w)  + C (w)r  (oj)  = Y (oj)  -C  (w)  v (w)  -C  (w)v  (w), 
sc  6 s c i c 2 s 


C (ai)r  (cj  ) - C (u)r  (u) 

5 S 6 C 


Y (d)  - C (aj)v  (w)  + C (w)v  (<u) 


and  solving  for  C (w) , we  find 


• » • . • • ••  •• 

C(r‘  + r2)  = Y r - Y r - c (v  r - v r ) 

6 c s'  cs  cc  1 c s sc 


(35) 


C (v  r + v r ) 
2 s s c c 


If  y(t)  is  the  position  of  the  forward  support  and  -y(t) 
the  position  of  the  after  support,  we  have 


and 


v ( t ) = - g2.  y ( t ) , r ( t)  = i y ( t) 

/v  U 

V (u ) = - -2-  V (u) 
c d c 

v (w)  = - ^ y (w) 
s d s 


I 


(36) 


_ £ *>  _ 


r (co) 
c 


£ y (co) 
a s 


rs(u)) 


- d yc(a,) 


An  examination  of  eqs.  (36)  shows  that  -v  r = v r , and 

^ s s c c 

therefore,  the  last  term  on  the  right  hand  side  of  eq.  (35) 

is  identically  zero  for  all  frequencies.  Substituting 

eqs.  (36)  into  (35)  we  find 


C (w) 
s 


d (Y  y - Y y ) - C ( co ) u (y  2 + y 2; 

sJ  s cx  c i o c s 

co  (y  “ + y 2 ) 
c 1 s 


For  the  case  of  zero  frequency,  it  is  a simple  matter  to 

show  that  both  Y and  y will  go  linearly  to  zero  as  co  goes 
s - s 

to  zero.  However,  both  Y and  y approach  non-zero  limits 

c c 

which  are  equal  to  the  difference  between  their  initial 
and  final  steady-state  values: 


Y ( oj=0  ) = Y_  - Y 
c T o 


yc(w=0)  = yT  - yQ 


Therefore,  when  co  = 0,  we  have 


lin  i Y - Y 

C (co=0)  = (o-o  — (-d  — _ c (co=0)u  + 0 (co2 ) ) 

6 w yT  " yo  1 ° 


Returning  to  the  equations  of  motion  (12),  we  can  see  that, 

uo 

for  the  steady-state  case  v = - 3—  y 
1 o d o 


and  r 


0, 
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b 


,oo 

V +v  / N ( x ) dx  = y 
yy  o o J o yy  o 


or 


(37) 

and  similarly 


u 

- C (w=0)  y 
i d o 


= Y 


:38) 


- c (u)=0 ) -H  v 

l d T 


= Y_ 


Finally  we  see  that  C (tu=0)  = 0.  However,  in  our  calcula- 

6 

tions  Ym/  Y , y . and  v are  all  measured  quantities  and 
it  is  apparent  that  an  error,  no  matter  how  small,  in  any 
of  these  quantities  will  cause  the  singular  behavior  ob- 
served in  the  stability  coefficients. 

Since  it  is  impossible  to  obtain  measurements  of  infin- 
ite accuracy,  the  following  scheme  was  adopted,  referred  to 
as  the  zero-frequency  correction.  Let  the  measured  value 

of  the  force  be  designated  by  Y and  the  value  of  the  noise 

1 m 

Y . Then  we  have 
n 

Y = Y + Y 
m n 

During  the  oeriod  prior  to  the  impulse  we  can  measure  Y 

m 

and  the  displacement  y (we  assume  that  we  have  the  capa- 
bility to  measure  the  displacement  with  greater  accuracy 
than  the  force).  Equation  (37)  gives  us  the  value  of  Y 

o 

for  this  period,  so  that  we  can  calculate  Y and  subtract 

n 

it  from  Y at  every  point  in  time.  We  now  make  a finite 
m 


« 
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Fourier  transform  of  Y (t)  - Y from  t=0  until  t=T , 

m n 

where  T is  the  time  when  Y (t)  - Y has  stabilized  at  a 

m n 

value  approximately  equal  to  Y^  as  defined  by  eq.  (38). 
The  assumption  is  made  that  Y is  exactly  equal  to  Y^  for 
all  time  t > T. 

In  terms  of  the  extended  definition  of  the  Fourier 
transform 


Y(u>)  = /T  Y (t)  e lti>tdt  - i 

J o w 


(V_i“T  - Yo> 


this  assumption  is  equivalent  to  the  replacement  of  the 
measured  values  of  the  force  before  and  after  the  impulse 
(which  necessarily  contain  some  error)  by  their  values  as 
calculated  by  eqs.  (37)  and  (38). 

It  should  be  noted  that  a similar  assumption  has  already 
been  made  for  the  case  of  pure  sway.  In  this  case  the 
assumption  is  that  the  force  must  be  identically  zero  both 
before  and  after  the  impulse. 

The  results  of  the  combined  sway  and  yaw  runs  were 
calculated  a second  time  using  this  zero-frequency  correction 
and  are  presented  in  Figures  29  through  36.  The  change  in 
the  damping  coefficients  is  slight  and,  finally,  we  obtain 
good  results  for  all  coefficients  in  the  zero- frequency  limit. 


A Test  on  the  Linearity  of  the  System 

The  linearization  scheme  which  led  to  the  equations  of 
motion  (12)  requires  that  the  lateral  and  angular  velocities 
be  "small"  in  comparison  with  the  forward  velocity.  If  this 


requirement  is  met,  the  experimental  results  should  be  ince- 
'fc  oendent  of  the  exact  nature  of  the  impulse  given  to  the  me  del. 

Since  both  the  lateral  and  the  angular  velocities  depend 
upon  the  peak-to-peak  amplitude  of  the  pulse,  a series  of 
experiments  was  run  using  various  amplitudes.  Manning  (1976) 
oresents  the  results  of  the  entire  series  of  experiments, 

with  amplitudes  of  0.40  to  1.00  inches  and  Froude  numbers  of 

I 

0.20  and  0.30,  using  both  the  full-pulse  and  the  step-pulse. 
Up  to  one  inch  amolitude  (the  maximum  oossible  with  our  PMM) 
no  systematic  variation  of  the  results  could  be  observed. 

It  would  appear  then,  that  the  one-inch  amplitude  does  not 
violate  the  linearization  assumption,  and  since  it  produces 
j the  best  signal- to-noise  ratio,  it  is  the  "preferred"  pulse. 

Figures  37  and  38  are  typical  results  of  this  test. 
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The  Prediction  of  Ship  Maneuvers 

Once  one  has  a complete  set  of  the  stability  coef- 
ficients for  a given  ship,  it  is  possible  to  predict  the 
lateral  and  angular  motions  of  the  ship  for  a given  set  of 
external  forces  and  moments.  Alternatively , one  could  pre- 
dict the  forces  and  moments  necessary  to  produce  a given 
path. 

In  the  present  study,  the  forces  and  moments  produced 
by  the  rudder  were  the  only  ones  considered.  The  problem 
then  becomes  one  of  finding  the  path  of  the  ship  for  a given 
rudder  command. 

The  Rudder 

Since  the  model  which  was  used  is  not  fitted  with  either 
a rudder  or  a propeller,  it  was  necessary  to  make  some 
assumptions  about  the  rudder  forces.  The  presence  of  a 
rudder  has  quite  a significant  effect  upon  the  overall 
stability  of  the  ship  and  the  selection  of  a particular 
rudder  can  cause  a radically  different  behavior  of  the  ship 
if  the  rudder's  contribution  to  the  damping  coefficients 
causes  the  ship  to  become  stable  rather  than  unstable. 
However,  if  we  consider  two  rudders,  both  of  which  lead  to 
a stable  ship,  the  predicted  maneuvers  will  differ  in  ab- 
solute value  but  not  in  their  general  behavior.  Therefore, 
if  we  are  careful  to  select  a rudder  which  is  large  enough 
to  insure  the  stability  of  the  ship,  we  will  be  able  to 


predict  the  general  behavior  of  the  ship  and  to  compare 
predictions  made  by  the  traditional  approach  with  the 
transient-motion  approach,  even  though  the  absolute  value 
of  the  predictions  may  differ  somewhat  from  predictions 
made  for  the  rudder  which  is  actually  on  the  ship. 

Therefore,  we  shall  assume  that  the  ship  is  outfitted 
with  a spade  rudder  which  has  an  aspect  ratio  of  two  and 
an  area  equal  to  2.2  per  cent  of  the  length  times  the  draft. 
Assuming  a taper  ratio  of  0.45,  the  dimensions  of  the  rudder 
(the  length  of  the  full-scale  ship  is  314.5  feet)  are: 


Area 

=' 

83  ft2 

Span 

= 

12.9  ft 

Section 

= 

NACA  0015 

Sweep  Angle 

= 

0 

Max.  Chord 

= 

8.9  ft 

Min.  Chord 

= 

4.0  ft 

where  the  maximum  chord  is  measured  at  the  intersection  of 
the  rudder  and  the  hull,  and  the  minimum  chord  is  measured 
at  the  tip  of  the  rudder. 

By  following  the  technique  recommended  by  Taplin  (1960) 
and  using  the  data  compiled  by  Whicker  and  Fehlner  (1953), 
it  is  found  that  for  small  rudder  anqles,  the  lateral  force 
exerted  on  the  hull  by  this  rudder  can  be  approximated  by 

Y„  = 3000  lb.  per  deqree  rudder  anale. 

E 

Taking  the  distance  between  the  center  of  gravity  and  the 
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rudder's  center  of  effort  to  be  144  feet,  we  approximate  the 
moment  by 

NE  = - 432,000  ft-lb  per  degree  . 

Change  in  Stability  Coefficients  Due  to  Rudder 

Since  our  experiments  were  run  using  a model  without 
either  a rudder  or  a propeller,  and  since  the  rudder  can 
contribute  significantly  to  the  added  mass  and  damping 
coefficients  of  the  ship,  it  was  necessary  to  add  a correc- 
tion term  to  the  experimentally  determined  coef f icients . 
Mandel  (1967)  suggests  a method  of  finding  the  correction 
terms.  Mandel  assumes  that  the  correction  terms  are  not 
functions  of  frequency  and,  of  course  we  would  rather  not 
make  this  assumption  since  we  wish  to  compare  predictions 
made  with  frequency-independent  coefficients  to  predictions 
made  with  frequency-dependent  coefficients.  However,  we 
would  expect  the  added  mass  and  damping  coefficients  of  a 
deeply  submerged  body  to  be  frequency-independent  and  there- 
fore the  frequency  dependence  of  the  rudder  correction  terms 
is  not  likely  to  be  too  great.  Therefore,  we  have  followed 
the  method  outlined  by  Mandel  and  have  reached  the  following 
results : 

AY'  = .0019 

v 

AY*'  = - . 00014 

v 

AN'  = + .00088 

v 


+ .000066 


A Y_ 


A Y ■ 


A N 


A N- 


+ . 00083 
+ .000066 

- .00040 

- .000030 


These  correction  terms  have  been  added  to  the  experi- 
mental results  and  the  (dimensional)  stability  coefficients 
for  the  full-scale  ship  were  calculated. 


Methods  of  Prediction 

One  of  the  standard  techniques  used  to  find  the  solu- 
tions of  eqs.  (12) , i.e.  the  equations  of  motion  which  allow 
for  "memory  effects",  is  the  Fourier  transformation.  The 
equations  of  motion  are  transformed  into  the  frequency  do- 
main, and  v((j)  and  r(ui)  , the  Fourier  transform^  of  v(t)  and 
r(t) , are  found  by  a f requency-by- frequency  solution  of  eqs. 
(15).  The  inverse  Fourier  transform  then  provides  us  with 
the  ability  to  find  v(t)  and  r(t). 

There  are  two  major  sources  of  error  in  such  a solution. 
The  first  source  is  simply  the  inaccuracy  inherent  in  the 
inverse  Fourier  transformation  of  a discrete  function  v(w). 
One  must  be  careful,  therefore,  to  choose  the  distance  Aw 
between  the  discrete  values  of  v(w),  to  be  sufficiently 
small.  The  second  source  of  error  is  more  difficult  to  con- 
trol. As  mentioned  earlier,  we  can  evaluate  the  stability 

coefficients  for  0 < u>  < to  where  w is  finite.  Therefore, 

— ~ i i 

infinite 


i 

mu 


we  must  replace  the 


integral  in  the  Inversion 


Theorem  by  a finite  integral  and  v(w)  must  be  of  such  a form 
that 

r00  ~ , , iwt , • -a.,  . . iuitj 
j v ( uj  ) e dw  - < > v (<•>)  e dm 
1 o 

Since  the  velocities  of  a snip  are  unlikely  to  co'tam  signi 
ficant  components  at  high  frequencies,  this  condition  was 
assumed  to  hold  with  sufficient  accuracy. 

The  traditional  equations  of  motion,  eqs.  (5),  are  much 
simpler  to  solve.  In  fact,  if  one  assumes  constant  coeffi- 
cients (stability  derivatives  evaluated  at.  o=0)  it  is 
possible  to  find  the  exact  solution  to  the  problem.  Further 
more,  it  is  also  possible  to  solve  the  equations  using  the 
Fourier  transformation.  Therefore,  if  we  solve  eqs.  (5) 
by  both  methods,  exact  and  Fourier  transform  solutions,  we 
shall  have  a measure  of  the  accuracy  of  the  computer  program 
that  calculates  the  inverse  F.urier  transform. 

A Check  on  the  Accuracy  of  the  Computations 

A computer  program  has  b>‘On  written  that  is  .'apable  of 
solving  both  the  traditional  equations  of  motion,  egs.  (3), 
(where  the  stability  derivatives  are  assumed  to  be  constant 
and  equal  to  their  zero-frequency  value)  and  the  transient- 
motion  equations  (12).  The  program  uses  the  Fourier  trans- 
form to  solve  both  sets  of  equations. 

A comparison  of  the  solutions  of  eqs.  (5) , as  computed 
by  the  program,  to  the  exact  solutions  of  eqs.  (5)  provides 
us  with  a check  on  the  accuracy  of  the  program  itself. 


■ 


Such  comparisons  were  made  for  two  different  rudder  commands. 
Letting  ' be  the  rudder  angle  in  degrees,  the  first  command 
was 

(0  t < 0 

<5  (t)  = < 1.5  t,  0 < t < 10 
1 [15  , t > 10 

i.e.  the  rudder  angle  is  increased  linearly  to  a maximum 
angle  of  15  degrees  in  ten  seconds.  The  second  rudder 
command  was  an  instantaneous  increase  in  rudder  angle  to 
the  same  15  degree  maximum. 


For  both  rudder  commands,  it  was  found  that  the  error 
in  the  computed  solution,  relative  to  the  exact  solution, 
was  less  than  two  per  cent.  In  terms  of  the  overall  accuracy 
of  the  experimentally  determined  stability  derivatives,  the 
accuracy  of  the  computer  program  is  quite  good. 


The  Predicted  Maneuvers 

The  two  rudder  commands  already  defined,  6 and  5 , 

1 z 

correspond  to  the  maneuver  known  as  the  turning  circle.  In 
addition  to  these  two  commands,  predictions  were  made  for  a 
simple  change  of  course,  5 and  for  the  initial  phases  of 

3 

a zig-zag  maneuver,  6 . Where 
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' 


T 


and 


5 (t) 


I 
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1. 5t  , 
15 

‘’5-1. 5 1 , 
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5 (t) 


0 

1 . 5 1 , 

15 
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0 
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20 
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60 
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< 

70 

t 

> 

70 

In  order  to  examine  the  effect  of  the  frequency  depen- 
dence, predictions  were  made  for  each  rudder  command,  based 
upon  the  traditional  equations  of  motion  (5)  as  well  as  the 
transient-motion  equations  of  motion  (12) . Let  v (t)  be  a 

5 

computed  solution  of  eqs.  (5) , v (t)  a comouted  solution  of 

\ 2 

eqs.  (12) , and  v (t)  an  exact  solution  of  eqs.  (5) ; then  we 
round  the  remarkable  result: 


| v (t)  - v ( t ) | < ! v ( t ) - v ( t ) 

5 12  5 e 


for  au  t and  for  each  rudder  command.  In  other  words,  the 
differences  between  predictions  based  upon  the  different 
equations  of  motion  were  always  less  than  the  error  inherent 
in  the  computational  technique,  which  is  quite  small. 

The  predictions  for  the  various  maneuvers,  based  upon 
eqs.  (12)  are  shown  in  Figures  39  through  42  . Pre- 

dictions based  upon  eqs.  (5)  arc  not  shown  since  the  differ- 
ences are  too  small  to  be  seen  graphically. 
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Conclusions 

Even  if  we  set  aside  the  question  as  to  which  set  of 
linearized  equations  of  motion  is  the  correct  one  and  under 
what  circumstances  it  is  correct  to  use  the  simpler  tra- 
ditional approach,  we  are  of  the  opinion  that  the  impulse- 
test  method  is  a superior  means  of  evaluating  the  stability 
derivatives . 

First  of  all,  the  impulse  technique,  even  if  repeated 
runs  are  made  to  avoid  the  previously  mentioned  difficulties, 
requires  substantially  fewer  experiments  than  the  regular- 
motion  technique  to  obtain  the  same  results. 

Secondly,  the  regular-motion  technique  is  incapable  of 
evaluating  the  stability  derivatives  at  very  low  frequencies, 
and  since  this  same  frequency  range  is  of  paramount  importance, 
static-sway  tests  and  rotating-arm  tests  (if  one  has  the 
facilities)  are  used  to  supplement  the  data.  However,  even 
if  these  additional  tests  are  performed,  they  assist  only  in 
the  evaluation  of  damping  coefficients  and  contribute  little 
to  the  evaluations  of  added  masses.  The  impulse  test  has 
the  capability  of  accurate  measurement  all  the  way  down  to 
zero  frequency. 

Third,  since  the  impulse  test  requires  measurement  of  data 
for  only  a very  limited  time,  it  becomes  a more  versatile 
method.  Problems  such  as  wave  reflection  can  be  avoided 
simply  by  recording  the  data  before  the  reflected  wave  has 
time  to  encounter  the  ship.  Similarly,  a short  towing  tank 
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ceases  to  present  a problem. 

However,  once  a complete  set  of  the  added  masses  and 
» damping  coefficients  has  been  determined,  it  would  appear 

that  the  choice  between  the  two  sets  of  linearized  equations 
of  motion  is  of  little  significance  in  the  prediction  of  ship 
maneuvers.  Even  the  idealized  case  of  instantaneous  rudder 
response,  a severe  test  of  the  effect  of  frequency  dependence, 
led  to  predictions  which  differed  only  slightly. 

Recent  papers  by  Fujino  and  Motora  (1975) , Nomoto  (1975), 
and  Fujino  (1975)  have  all  expressed  the  opinion  that  the 
memory  effect  is  often  small  and  that  the  traditional  equa- 
tions of  motion  will  yield  good  results  for  many  standard 
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jure  9:  Results  of  regular-motion  tests  and  full-pulse  tests  for  Fn 


Figure  10:  Results  of  regular-motion  tests  and  full-pulse  tests 
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Figure  25:  Results  from  5 step-pulse  runs. 
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Results  from  5 step-pulse  runs,  with  zero-correction. 


Figure  36:  Zero-corrected  results  from  5 step-pulse  runs. 
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